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Abstract 

In this paper, we investigate the liquid gas phase transition for the spin polarized nuclear matter. 
Applying the lowest order constrained variational (LOCV) method, and using two microscopic 
potentials, AVis and C/R14+TNI, we calculate the free energy, equation of state, order parameter, 
entropy, heat capacity and compressibility to derive the critical properties of spin polarized nuclear 
matter. Our results indicate that for the spin polarized nuclear matter, the second order phase 
transition takes place at lower temperatures with respect to the unpolarized one. It is also shown 
that the critical temperature of our spin polarized nuclear matter with a specific value of spin 
polarization parameter is in good agreement with the experimental result. 
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INTRODUCTION 


Through the formation of protoneutron stars and in the relativistic heavy ion collisions, 
the symmetric nuclear matter can reach the high temperatures about 20 — 50 MeV l|, l2|. 
Therefore, liquid gas phase transition in nuclear matter may occur. To investigate the na¬ 
ture of liquid gas phase transition and the properties of critical point, many experimental 




ing a phenomeno 
about 17.5 MeV 


y-Gl and theoretical |7Hl4| studies have been done. The extracted value of critical temper- 

n 

ature from the heavy ion collisions is Tc ~ 16.6 MeV [ 6 (. However, the critical temperature 
driven by theoretical methods depends on the applied models and interactions. Using the 
variational method based on the UVu two-body nucleon-nucleon (NN) interaction includ- 
ogical three-body force (TNI), the critical temperature has been gotten 
. In the frame of relativistic Brueckner-Hartree-Fock (BHF) theory and 
applying the potentials given by Brockmann and Machleidt, it has been shown that the 
critical temperatures are smaller than those of the nonrelativistic investigations 8 |. Ap¬ 
plying different versions of scalar derivate coupling suggested by Zimanyi and Moszkowski, 
the critical temperature has been obtained in the range of 13.6 — 18.3 MeV 9|. Using 
the Bloch-De Dominicis diagrammatic expansion and employing the AU 14 as the bare NN 
interaction, without and with a phenomenological three-body force, the gotten critical tem¬ 
peratures are — 21 MeV and Tc ~ 20 MeV, respectively 1^. In the framework of hnite 
temperature BHF approach extended to include the contribution of a microscopic three- 
body force, it has been shown that the three-body force reduces the critical temperature 


from Tc ~ 16 MeV to 13 MeV 


11| . It has been seen that the critical temperature for the 


CDBONN potential in the BHF approximation {Tc — 23.3 MeV) is larger than the critical 
temperature for the same potential but in the Self-Consistent Greens Functions approach 
(SCGF) (Tc ~ 18.5 MeV) 12|. In this reference, it has also been shown that for AVig 
potential, the critical temperature is lower than for the CDBONN potential 12|]. Employ¬ 
ing the self-consistent Hartree-Fock approach using different mean-held interactions of the 
Skyrme and the Gogny types, it has been indicated that the critical temperatures span a 
wide range of values, from Tc ~ 14 to 23 MeV showing the effective interaction dependence 
of the critical properties . In addition, the predicted value of critical temperature for the 


2 




















density-dependent relativistic mean-field models is abont 15.7 MeV [l^ . 

Previonsly, we have investigated the liqnid gas phase transition of spin nnpolarized asym¬ 
metrical nnclear matter, nsing the lowest order constrained variational (LOCV) method 15|. 
The dependence of the critical qnantities on the isospin polarization of nnclear matter has 
been discnssed. In the present work, we are interested to consider the effect of spin polariz¬ 
ability of symmetric nnclear matter on the liqnid gas phase transition applying the LOCV 
calcnlation nsing the microscopic potentials. 


II. FINITE TEMPERATURE CALCULATIONS FOR SPIN POLARIZED NU¬ 
CLEAR MATTER WITH THE LOCV METHOD 


We consider a system of A interacting nncleons. This system is composed of spin-np and 
spin-down nncleons with densities p^\ p^'> and Labels (-I-) and (—) are nsed for 

the spin-np and spin-down nentrons (n) and protons (p), respectively. The total densities 
for nentrons (pn), protons (pp) are given by, 

o = -I- 

Fn Fn ' Fn 1 

Pp = ■ ( 1 ) 

The total density of system is 


P — Pn + Pp- 


( 2 ) 


In the case of symmetric nnclear matter, we have pn = Pp. The spin polarization parameters 
which describe the spin asymmetry of the system are defined as. 


6p = 


Pi~^^ - pi ^ 


= 


p{+) — pi-) 
rn rn 


( 3 ) 


Pp pn 

To obtain the macroscopic properties of this system, we shonld calcnlate the total free 
energy per nncleon, F, 


F = B - T{SiP + S'pl + Sip + S<->), 

(i) 


( 4 ) 


where E is the total energy per nncleon and Sj ' is the entropy per nncleon corresponding 
to the isospin and spin projection j, i, respectively. 


Sfp.T) = -li:{[l-nf(i,r,p<‘')]lii[l-nf(i,r,pf) 


A 


+r^ {k,T, pf) In nf {k,T, pf )} . 


( 5 ) 
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In above equation, m \k,T, Pj) is the Fermi-Dirac distribution function, 






( 6 ) 


exp (^[ef{k,r,pf) - pf{T,pf)] /ksT) + l’ 

where is the single particle energy of a nucleon, and is the chemical potential which is 
determined at any adopted value of the temperature T and number density by applying 
the following constraint, 

i:nf(fc,r,pf) = ivf. (7) 

k 

In our formalism, the single particle energy of nucleons with momentum k approximately is 
written in terms of the effective mass as follows 


e*^(fc,r,p*0 = 


'k) 


+ f/«(r,p®). 


( 8 ) 


where p*-*^) is the momentum independent single particle potential. In fact, we use a 

quadratic approximation for the single particle potential incorporated in the single particle 
energy as a momentum independent effective mass. We introduce the effective masses, m 


as variational parameters 


A 


ii) 
j ’ 

18l | . We minimize the free energy with respect to the variations 


in the effective masses and then we obtain the chemical potentials and the effective masses 
of the spin-up and spin-down nucleons at the minimum point of the free energy. This 
minimization is done numerically. 

For calculating the total energy of polarized symmetric nuclear matter, we use the LOCV 
method. We choose a trial many-body wave function of the form 


'ip = 


( 9 ) 


where 0 is the uncorrelated ground state wave function (the Slater determinant of plane 
waves) of A independent nucleons and = J^(l ■ ■ • A) is an appropriate A-body correlation 
operator which can be replaced by a Jastrow form i.e.. 


^ = 511 /®). 

i>j 


( 10 ) 


where iS is a symmetrizing operator. Now, we consider the cluster expansion of energy 
functional up to the two-body term |i^ . 


Pdf]) _ 1 i'^W) -PIP 


( 11 ) 
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For the hot nuclear matter, the one-body term Ei is 


27.2 


^1= i: i: sM-. 

j=n,pi=+- k 


nf{k,T,pf). 


( 12 ) 


The two-body energy E 2 is 


where 


^2 = 1^(12)hi-i*), 

ij 


>'(12) = -|(^l/(12). [V;2./(12)1] + f(12)V(12)f(12). 


(13) 


(14) 


In above equation, /(12) and l'^(12) are the two-body correlation function and potential. In 


our calculations, we use the microscopic potentials {AVi^ |20(] and 171h4+TNI |22|). The two- 


body correlation function, /(12), which is induced by the strong force is given by /(12) = 
Yl\=i where Pi^ has been given in Ref. 21|. Using this two-body correlation 

function and the microscopic potentials, after doing some algebra, we get an equation for 
the two-body energy. In the next step, we minimize the two-body energy with respect to 
the variations in the functions subject to the normalization constraint, ~ 

P{12)\ij)a = 0, 2l|]. In the case of polarized symmetric nuclear matter, the Pauli function 
hs^{r) is as follows [21 1 


hsAr) = 


1 


-1/2 


; 5'^ = ±1 
; >^. = 0 


(15) 


where 




' n‘'‘>(k,T,p<‘^)Jo(kr)k^-dk. 


(16) 


From the minimization of the two-body cluster energy, we get a set of coupled and uncoupled 
differential equations 0. By solving these diffeientia. e,nat.ons. we ean obtain coiielation 
functions to compute the two-body energy term. For more details see Refs j^, 


III. RESULTS AND DISCUSSION 


In this paper, we study the liquid gas phase transition for the nuclear matter in the 
unpolarized, ferromagnetic and antiferromagnetic states. It should be mentioned that the 
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cases Sp = 6n = 0, 6p = dn and 6p = —6r,. are called the unpolarized, ferromagnetic and 


antiferromagnetic states, respectively 
may also called polarized states. 


— 

5n 8 

23, 

241 


24l |. Ferromagnetic and antiferromagnetic states 



FIG. 1: The free energy per particle versus the density with [/F 14 +TNI and AVig potentials for 
different magnetic states (with 6p = ±0.5 for the ferromagnetic state and also 6p = ±0.5 for the 
antiferromagnetic state) at a fixed temperature. 

Fig. [1] shows the free energy per particle versus the density with f/hi 4 ±TNI and AVig 
potentials for different magnetic states. Obviously, at lower densities, the convexity of free 
energy which is the condition of stability, is violated and nuclear matter is mechanically 
instable. This leads to a first order liquid gas phase transition in nuclear matter. From Fig. 
[H it can be seen that at higher densities, the convexity condition remains valid, and the 
system is stable. We can see that in the case of t/Vi 4 ±TNI potential, the free energy is more 
sensitive to the variation of density. 

Fig. |2]presents the pressure-density isotherms at a hxed temperature with 17Fi4±TNI and 
AVis potentials for three magnetic states. We see that for all magnetic states, the isotherms 
corresponding to temperature T = 13 MeV are below those of the corresponding critical 
temperatures. At different magnetic states, the system shows a mechanical instability. For 
both applied potentials, in the antiferromagnetic state, we see the stiffest equation of state 
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FIG. 2; Same as Fig. [T]but for the pressure. 

(EOS) for a fixed temperature. This leads to the lower value of the critical temperature 
and critical density. By comparing the isotherms, we realize that for unpolarized state, the 
instable region is the most extensive. This indicates that for unpolarized system, the second 
order phase transition occurs at a higher temperature. 

Figs. |3] and H] show the pressure-density isotherms at different temperatures with 
t/I/i 4 -|-TNI and AVig potentials for three magnetic states. The isotherms present a typi¬ 
cal Van der Waals like behavior in which the liquid and gaseous phases coexist. For all 
magnetic states, we have found that the instable region reduces by increasing the tem¬ 
perature. The properties of the mixed phase in thermal equilibrium can be calculated by 
applying the equal-area Maxwell construction. To calculate the pressure of nuclear matter 
during the phase transition, we employ the Maxwell construction as shown in Figs. |5]and 
El Between the left and right crossing points, the nuclear matter does not experience the 
unstable curve (the solid line), but physically it goes through the curve at constant pressure 
(the dashed line). The stable states of nuclear matter between the left and right crossing 
points correspond to the points on the dashed line. In fact, the continuous evolution with 
phase mixing appears instead of the discontinuity of the EOS known as the hrst order phase 
transition. 
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Fig. [7] shows the coexistence curves, the liquid and gas densities versus temperature, for 
t/Vi 4 +TNI and AVis potentials at different magnetic states. These curves are determined 
by the Maxwell construction. As the temperature increases, the liquid density decreases 
and the gas density increases. Inside the coexistence curve, the stable phase of nuclear 
matter is a mixture of liquid and gas. The coexistence region decreases by increasing the 
temperature. This region disappears at the critical temperature, and the densities become 
equal in which the nuclear matter experience a second order phase transition. For polarized 
states at each temperature, the coexistence region is smaller than in the case of unpolarized 
state. In addition, for the antiferromagnetic state, the coexistence region is more limited 
than the ferromagnetic one. With two employed potentials, the coexistence region for the 
antiferromagnetic state disappears before the ferromagnetic one. This indicates that the 
critical temperature is lower for the antiferromagnetic state. At each temperature by com¬ 
paring the coexistence regions obtained using two employed potentials, it is obvious that 
the coexistence region is wider for AVis potential. At the critical temperature, Tc, with the 
condition {^)t, 5 p = {^^)t, 5 p = 0 , the coexistence region changes to a point corresponding 
to the critical pressure Pc and critical density pc- Our results for the critical temperature 
are presented in Fig. [ 8 l We can see that the critical temperature is a symmetric function of 
the spin polarization parameter. This quantity decreases by increasing the magnitude of the 
spin polarization parameter. The effect of spin polarization parameter on the critical tem¬ 
perature is more signihcant for AVis potential. Moreover, the effect of the spin polarization 
parameter on the critical temperature is more signihcant for the antiferromagnetic state. At 
each value of the spin polarization parameter, the critical temperature for the ferromagnetic 
state is greater than the antiferromagnetic one. The difference in the critical temperatures 
of the ferromagnetic and antiferromagnetic states is signihcant for larger values of polariza¬ 
tion. Using Fig. [HI it is possible to extract the value of the spin polarization parameter at 
which the corresponding critical temperature is in agreement with the experimental value, 
Tc — 16.6 MeV. For the ferromagnetic nuclear matter with [/U 14 -I-TNI and AVis poten¬ 
tials, the critical temperatures are equal to 16.6 MeV with 6p = ±0.31 and 6p = ±0.44, 
respectively. In addition, for the antiferromagnetic nuclear matter with 17Vi4±TNI and 
AVis potentials, the second order phase transition occurs at 16.6 MeV when 5p = ±0.25 
and 5p = ±0.41, respectively. 

Our results for the critical density and pressure are presented in Tables HI and [TTl Table |T] 


indicates that the critical density for the case 6p = ±0.25 is larger than that of 6p = ±0.50. 
In addition, with AVi^ potential, the critical density is more sensitive to the spin polarization 
parameter. It should be noted that a similar behavior can be seen for the critical pressure 
from Table Ull 

In the case of liquid gas phase transition, the difference in the liquid and gas densities 
plays the role of order parameter. We dehne the order parameter as m = puquid — Pgas 
to study the critical properties of the nuclear matter. This quantity is presented in Fig. 
m for the 17 I/i 4±TNI and AVig potentials at different magnetic states. It is obvious that 
the order parameter vanishes at the critical temperature. At each temperature, below the 
critical point, for the 17Vi4±TNI potential, the order parameter at the polarized states is 
smaller than that of the unpolarized state. It is evident that the order parameter for the 
ferromagnetic state is larger than the antiferromagnetic one. 

Fig. [TT] presents the entropy per particle at the critical pressure as a function of tem¬ 
perature. We can see that the entropy continues at the phase transition point, while as we 
will see, the heat capacity diverges with a power law behavior at the phase transition point 
indicating that the phase transition is of second order. At each temperature, the entropy of 
nuclear matter in antiferromagnetic state is larger than in the ferromagnetic and unpolarized 
one. The heat capacity at the critical pressure of the nuclear matter, Cp = T(|^)p, is plotted 
in Fig. [121 This quantity diverges with a power law behavior at the critical temperature 
conhrming that the phase transition is the second order. 

Fig. [13] shows the isothermal compressibility, Kt = s-s a function of temper¬ 

ature. Below Tc, we have shown the liquid branch. For all temperatures, the isothermal 
compressibility is also computed at Pc- It is obvious that near the critical temperature from 
either side, the isothermal compressibility diverges with a power law behavior. This phase 
transition is specihed W an order parameter which is non zero below the critical temperature 
and is zero above it , as we saw in Fig. [TOl Close to the critical point, the fluctuations 
intervenes dominate |26l |. 

The critical temperature obtained using different potentials and methods have been com¬ 
pared in Table IIIII It can be seen that for AVis potential with SCGF method and also 
for AVis + TNI potential with BHF method, the extracted critical temperatures are lower 
than those of our results. Among the results given by UVu ± TA^J, AI /14 and AVu ± TA^J 
potentials with the variational calculations and BlochDe Dominicis expansion, it is obvious 
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TABLE I: Our results for the critical density of symmetric nuclear matter for ?7Vi4+TNI and AEig 
potentials at different magnetic states. 


Potential 

Magnetic State 

Sp 

Pc 

UVu+TNl 

Unpolarized 

0.00 

0.07 

Ferromagnetic 

±0.25 

0.07 

±0.50 

0.06 

Antiferromagnetic 

±0.25 

0.07 

±0.50 

0.06 

00 

Ferromagnetic 

±0.25 

0.13 

±0.50 

0.11 

Antiferromagnetic 

±0.25 

0.13 

±0.50 

0.09 


TABLE II: Same as Table |T] but for the critical pressure. 


Potential 

Magnetic State 

Sp 

Pc (MeVfm-^) 

uVu+Tm 

Unpolarized 

0.00 

0.28 

Ferromagnetic 

±0.25 

0.21 

±0.50 

0.13 

Antiferromagnetic 

±0.25 

0.19 

±0.50 

0.11 

00 

Ferromagnetic 

±0.25 

0.82 

±0.50 

0.30 

Antiferromagnetic 

±0.25 

0.76 

±0.50 

0.33 


that the value ~ 21 MeV corresponding to AV 14 and BlochDe Dominicis expansion is 
almost close to our result, ~ 24.3 MeV [1^. Our critical temperature for the unpolarized 
system, — 24.3 MeV, is in a nearly good agreement with the obtained valne by the 
CDBONN potential in the BHF approximation (Tc ~ 23.3 MeV). 
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TABLE III: Critical temperature for different potentials and methods (with 6p = ±0.5 for our 
polarized states). 


Potential 

Method 

Tc (MeV) 

Reference 

17Vi4±TNI 

LOCV 

18.0 

Our result(Unpolarized) 

17Vi4±TNI 

LOCV 

15.1 

Our result (Ferromagnetic) 

17Vi4±TNI 

LOCV 

14.1 

Our result (Antiferromagnetic) 

00 

LOCV 

15.0 

Our result (Ferromagnetic) 

00 

LOCV 

14.0 

Our result (Antiferromagnetic) 

AFi8 

LOCV 

24.3 

[la 

AFi8 

SCGF 

11.6 

[12] 

00 

BHF 

18.1 

[12] 

CDBONN 

SCGF 

18.5 

[12] 

CDBONN 

BHF 

23.3 

[12] 

17Vi4±TNI 

Variational Calculation 

17.5 

m 

AVi4 

BlochDe Dominicis expansion 

21 

[10] 

AVi4±TNI 

BlochDe Dominicis expansion 

20 

[10] 

AVi8±TNI 

BHF 

13 


Experimental value 


16.6 

m 


IV. CRITICAL EXPONENTS FOR SPIN POLARIZED NUCLEAR MATTER 

In this paper, we are interested in the critical exponents of spin polarized nuclear matter 
which describes the behavior of the thermodynamic properties of our system near the critical 
point. We consider the critical isotherm, 

F-Pc ~ (p-Pc)^;p —> Pc, (17) 

where T = T^. The slope of this thermodynamic property on log-log scale corresponds to the 
critical exponent S. The values of this critical exponent with FVi 4 ±TNI and AVis potentials 
for different magnetic states which are calculated using Fig. [9] are presented in Table HVl 
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By definition of 


_ T-T, 

^ T ’ 

and considering the order parameter as a function of £ on the log-log scale in Fig 
can calculate the exponent (3 for this parameter, 

m = ^ 0. 

Our results for the value of /3 with f/l/ 14 -l-TNI and AVis potentials for different magnetic 
states are presented in Table HVl 

The heat capacity is related to the exponent a by 

Cp = (—e) “;£—)■ 0 . ( 20 ) 

The values of the critical exponent a with f 7 Vi 4 -|-TNI and AVis potentials for different 
magnetic states which are calculated using Fig. [12] are presented in Table HV] 

To describe the isothermal compressibility near the critical point, we calculate the critical 
exponent 7 , 


(18) 
. EUl we 

(19) 


Kt = 0 . 


( 21 ) 


The values of this critical exponent which are obtained by Fig. [13] are also presented in 
Table IIVI From Table IIVI we can see that near the critical point, the critical exponents 
calculated with 17Vi4-|-TNI and AVis potentials for different magnetic states are nearly 
identical. However, there are some small differences in the values because of the numerical 
errors occurring in the calculations near the critical point. From Table [IV] it is evident that 
the Griffiths and Rushbrooke inequalities j^. 28| 


a + 2/3 + j > 2, 


( 22 ) 


and 


cx -h /3[1 -l- 5) ^2, 


(23) 


are established for our critical exponents. 


12 






TABLE IV: The critical exponents for the symmetric nuclear matter with [/V 14 +TNI and AVig 
potentials for different magnetic states (with 5p = ±0.5 for the polarized states). 


Critical exponent 

Potential 

Magnetic State 


d 

C/Vi4 ±TNI 

Unpolarized 

0.60078 ± 0.02533 

Ferromagnetic 

0.46184 ±0.01729 

Antiferromagnetic 

0.55648 ±0.01592 

ALis 

Ferromagnetic 

0.53167 ±0.01598 

Antiferromagnetic 

0.59037 ±0.00703 

(5 

C/Vi4 ±TNI 

Unpolarized 

2.54265 ± 0.04097 

Ferromagnetic 

2.65974 ±0.03192 

Antiferromagnetic 

2.59537 ± 0.02929 

AVis 

Ferromagnetic 

2.59047 ±0.03184 

Antiferromagnetic 

2.59117 ±0.04159 

a 

uVu+Tm 

Unpolarized 

0.68480 ± 0.05970 

Ferromagnetic 

0.66717 ±0.03961 

Antiferromagnetic 

0.60539 ±0.04601 

AVis 

Ferromagnetic 

0.76092 ± 0.00636 

Antiferromagnetic 

0.81497 ±0.01300 

1 

f7Vi4±TNI 

Unpolarized 

0.77925 ± 0.00334 

Ferromagnetic 

0.87677 ±0.01601 

Antiferromagnetic 

0.86563 ±0.02087 

AV 18 

Ferromagnetic 

0.81733 ±0.01216 

Antiferromagnetic 

0.80605 ±0.00806 


V. SUMMARY AND CONCLUSIONS 

In this work, we have calculated the thermodynamic properties of spin polarized nuclear 
matter to study the liquid gas phase transition using the lowest order constrained varia¬ 
tional method employing the microscopic potentials. The critical properties of unpolarized, 
ferromagnetic and antiferromagnetic nuclear matter have been considered. By applying the 
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Maxwell construction to study the mixed phase in thermal equilibrium, we have determined 
the pressure of nuclear matter during the phase transition and also the coexistence curve. We 
have seen that the extension of the coexistence region is different for three magnetic states. 
For unpolarized system, the second order phase transition occurs at a higher temperature. 
In addition, for each value of the spin polarization parameter, the critical temperature is 
greater for the ferromagnetic state. Studying the critical temperature for different spin 
polarization parameters, it has been shown that the critical temperature is a symmetric 
function of the spin polarization parameter. Moreover, an increase in the magnitude of the 
spin polarization parameter reduces the critical temperature. In investigation of the critical 
temperature, it has been found that the antiferromagnetic nuclear matter is more sensitive 
to the value of the spin polarization parameter. The critical temperature of our polarized 
nuclear matter with a specihc value of spin polarization parameter is in a good agreement 
with the experimental result. Studying the order parameter for the liquid gas phase tran¬ 
sition, we have seen that below the critical point, the order parameter for the unpolarized 
nuclear matter is larger than that of polarized one. It has been also clarihed that the order 
parameter for the ferromagnetic state is larger than the antiferromagnetic one. For our 
system, the heat capacity and isothermal compressibility diverge with a power law behavior 
at the phase transition point. This confirms that the phase transition is of second order. 
Finally, we have calculated the critical exponents of the spin polarized nuclear for different 
magnetic states, concluding that the Griffiths and Rushbrooke inequalities are established 
for these exponents. 
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FIG. 3: The pressure-density isotherms for [/F 14 -I-TNI potential at different temperatures for three 
magnetic states (with 6 p = ±0.5 for the polarized states). 
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FIG. 4: Same as Fig. [3] but for AVig potential. 
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FIG. 6: Same as Fig. [5] but for AVi% potential. 
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FIG. 7: The liquid gas coexistence curves for ?7 Fl4 +TNI and AVis potentials at different magnetic 
states (with Sp = ±0.5 for the polarized states). 
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FIG. 9; The {P — Pc) versus {p — pc) at critical temperature (Tc) on log-log scale with 17Ti4+TNI 
and AVis potentials for different magnetic states (with 6p = ±0.5 for the polarized states). 
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FIG. 10: Two left panels; The order parameter for the liquid gas phase transition as a function 
of temperature with 17Vi4+TNI and AVis potentials for different magnetic states (with Sp = ±0.5 
for the polarized states). Two right panels; The order parameter versus e on log-log scale. 
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FIG. 11: The entropy per particle as a function of temperature at critical pressure with t/Vi 4 +TNI 
and AVis potentials for different magnetic states (with 5p = ±0.5 for the polarized states). 
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FIG. 12: Two left panels: The heat capacity per particle at critical pressure as a function of 
temperature with t/Fi 4 +TNI and AVis potentials for different magnetic states (with 6p = ±0.5 for 
the polarized states). Two right panels: The heat capacity per particle at critical pressure versus 
6 on log-log scale. 
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FIG. 13: Two left panels: The isothermal compressibility as a function of temperature with 
t/Fi 4 +TNI and AVig potentials for different magnetic states (with 5p = ±0.5 for the polarized 
states). Two right panels: The isothermal compressibility versus e on log-log scale. 
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